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Summary. In this chapter we present a general method of constructing the effective equation which 
describes the behaviour of small-amplitude solutions for a nonlinear PDE in finite volume, provided 
that the linear part of the equation is a liamiltonian system with a pure imaginary discrete spectrum. 
The effective equation is obtained by retaining only the resonant terms of the nonlinearity (which may 
be hamiltonian, or may be not); the assertion that it describes the limiting behaviour of small-amplitude 
solutions is a rigorous mathematical theorem. In particular, the method applies to the three- and four- 
wave systems. We demonstrate that different possible types of energy transport are covered by this 
method, depending on whether the set of resonances splits into finite clusters (this happens, e.g. in case 
of the Charney-Hasegawa-Mima equation), or is connected (this happens, e.g. in the case of the NLS 
equation if the space-dimension is at least two). For equations of the first type the energy transition to 
high frequencies does not hold, while for equations of the second type it may take place. In the case 
of the NLS equation we use next some heuristic approximation from the arsenal of wave turbulence to 
show that under the iterated limit “the volume goes to infinity”, taken after the limit “the amplitude 
of oscillations goes to zero”, the energy spectrum of solutions for the effective equation is described by 
a Zakharov-type kinetic equation. Evoking the Zakharov ansatz we show that stationary in time and 
homogeneous in space solutions for the latter equation have a power law form. Our method applies to 
various weakly nonlinear wave systems, appearing in plasma, meteorology and oceanology. 


1.1 Introduction 

It is well known that solutions of linear evolution PDEs in finite volume are superpositions of normal 
modes of oscillations (in most cases of interest these are the Fourier modes). When a nonlinearity is 
added as a perturbation, different modes start to interact and the solutions of the equation can be 
approximated by suitable power series expansions, provided that the nonlinearity is sufficiently small 
(or, in other words, the PDE is weakly nonlinear). In such cases, the equation can be written as 

ut + L(u) = eN(u) , (1.1) 

where L is the linear operator, N denotes the nonlinearity and e is a small parameter, 0 < e 1. The 
equation may contain a stochastic force, and in that case it reads 


Ut + L(u) = eN(u) + \fs (random force) 


( 1 . 2 ) 
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(the scaling of the random force by the factor \fe is the most natural, see below). We will show that 
the limiting, as e —> 0, exchange of energy between the modes may be described by replacing the 
original system with a suitable effective equation. This result may be regarded as a PDE-version of 
the Bogolyubov averaging principle (see [1]) which implies a similar property for distribution of energy 
between the oscillating modes for small-amplitude oscillations in finite-dimensional nonlinear systems. 

The mentioned above convergence that holds as e —» 0 and various properties of the corresponding 
effective equations have been rigorously established (see [2, 3, 4, 5, 6, 7] and the discussion in [7]). 
The treatments of the deterministic and stochastic equations are similar, but the results, obtained in 
the stochastic case, are significantly stronger: while the deterministic effective equation controls the 
dynamics only on time intervals of order e _1 , in the presence of stochastic forcing the corresponding 
effective equation also approximates the stationary measure for the original equation, thus controlling 
the asymptotical in time behaviour of solutions when e <C 1. Moreover, in the absence of forcing we only 
get information concerning the exchange of energy between the modes, whereas in the stochastically 
forced case, the stationary measure for the effective equation controls both the energies and the phases 
of the normal modes of solutions. 

Below we explain how to construct the effective equations for eq. (1.1) and eq. (1.2) from the 
resonant terms of the nonlinearities. We will discuss two examples: the nonlinear Schrodinger equation 
and the Charney-Hasegawa-Mima equation on the f) plane. These two equations display completely 
different types of energy exchange between modes, and we will explain why this happens. 


1.2 How to construct the effective equation 

We consider hamiltonian PDEs, whose linear parts have imaginary pure point spectra and are diagonal 
in Fourier modes. Written in terms of the complex Fourier coefficients v = {u k } (also called waves), 
the equations which we study read 

■^Vk = iuj k v k + Pk{v) , k £ Z d . (1.3) 

Here oj k are real numbers and P(v) = (Pk(v), k £ Z d ) is a polynomial nonlinearity in v of certain order 
q , of the form 

Pk(v) = J2 Ck 1 ... kgk Vk 1 • • • « kj Xp+i ''' w k 9 <5p+i P .., k , 

p<q k i... k p kp+1 ... kg 

where c kl ... k<jk are some complex coefficients, v* is the complex conjugate of v and 

ci...p _ J 1 if ki + ... + k p = kp + i + ... + kq + k . . 

p+ i~« k - \ 0 else ' 1 ' ’ 

We always assume that “the nonlinearity does not pump energy in the system”: 

Re ^2 Pk(v)vk < 0 (1.5) 

k 

(in most case of interest the l.h.s. vanishes). 

The quantities J k := |u k | 2 /2, Ek = w k / k and <fik = Arg(v k ) are called, respectively, the wave action, 
wave energy and wave phase. The relation between u and k, i.e. the function k —> w k , is called the 
dispersion relation, or dispersion function. 
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The weakly nonlinear regime corresponds to solutions of small amplitude e. We will study it in 
the presence of damping and, possibly, a random force, whose magnitude is controlled by another 
parameter, call it v. So, instead of (1.3), we will consider 

= «u k «k + e q P k (v) - I'TkVk + irVnb k f3 k , k G TL d , (1.6) 

where 7k > 7» > 0 controls the damping term, 6k > 0 controls the forcing and the parameter /r £ {0,1} 
is introduced to consider at the same time the forced and non-forced cases. The /3 k ’s are complex white 
noises, independent from each other. 3 The factors v and yjv in front of the damping and the dissipation 
are so chosen that, in the limit of v —> 0, the solutions stay of order one, uniformly in t > 0. 

Note that, while £ controls the size of the solutions, 1/zz is the time-scale on which the forcing acts 
significantly, as it is the time needed for the standard deviations of the processes yT73 k to become of 
order one. If /i = 0 and the system (1.6) is deterministic, still its time-scale is l/v since, as we explain 
below, its solutions with initial data of order one stay of order one for t < z/ _1 , while for much bigger 
values of time they are very small since in view of (1.5) their ^ 2 -norms satisfy 

K*)ll 2 < K0)|? 2 e“‘ ! ' 7 * . 


We will consider the regime 


(where q is the degree of P), and study solutions of the equation with given initial conditions on the 
time-scale l/v, examining them under the limit v —> 0. Passing to the slow time r = vt (so that time 
i ~ l/v corresponds to r of order 1), eq. (1.6) becomes 

Uk = *zz _1 w k Uk + Pk{v) - 7 kVk + nbk/3 k , k e Z d , (1.8) 

where the upper dot stands for 4~. 

We claim that, in the limit when v (or, equivalently, e) goes to zero, the distribution of the energies 
E k on times r of order one is described by an effective equation whose nonlinearity is constituted by 
resonant terms of the nonlinearity (see (1.13) below). 

It is easier to understand the role of resonances and the form of the effective equation by passing to 
the interaction representation ( cf. [1, 8, 9]), i.e., by performing the time-dependent change of variables 
from «k to 

«k = e-^ _1 “ kT Uk , (1.9) 

which transforms (1.8) to 

h k = Rk(a, z^ _1 r) — 7k(ik + qib k e~ lI/ “ kT /3 k , k € Z d , (1.10) 


where Rk denotes the nonlinearity, written in the a-variables. That is 


R k (a,z/V)=^ c ki...k,kU kl •••Uk p Uk p+ i---Uk,<5p +1 p .. 9k 

p<q k l ...kpk p+1 ...k, 

x exp (ii/~ 1 r (tu kl + • • • + tukp — Wk p+1 — ... — tUk, — w k )) • 

Noting that the collection of the processes {e~ v “ kT /3k} is another set of standard independent complex 
white noises, we re-write eq. (1.10) as 

3 That is, /3 k = ( d/dt)f3 k , /3 k = f3 k +if 3k~, where /3 k are standard independent real Wiener processes. 
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d k = Rk(a, v r) - 7 k ak + /x&k/3 k , k G Z d . 
In the sum defining Rk, the terms for which the resonance conditions 


Wki + • • • + L0 kp — Wk p+1 + . . . + LO kq + CJk 
ki + • • • + k p — kp+i -(-••• kq + k 


( 1 . 11 ) 


( 1 . 12 ) 


are satisfied (called the resonant terms) under the limit v —» 0 behave completely differently from 
others terms (called the nonresonant terms). Namely, the nonresonant terms oscillate faster and faster, 
whereas the resonant terms do not. We will say that a set of Z d -vectors {ki,... ,kq,k} forms a reso¬ 
nance if relations (1.12) are satisfied, if Ck^.k^k ^ 0, and the set {ki,... , k p } does not equal the set 
{kp+i,..., kq, k}. The collection of all resonances is called the resonant set. 

In the spirit of the finite-dimensional averaging, following the Bogolyubov averaging principle (see 
[1], the behaviour of solutions of (1.11) under the limit v —> 0 is obtained by replacing the nonlinearity 
Rk with its time average, i.e. with 

1 f T 

Inn — / R k (a,t)dt . 

T—>oo T J q 

Since for any real number A we have 


lim 

T —>oo 


l 

T 


iXt 71 

e at 


1 if A = 0 
OifA^O ’ 


then only the resonant terms survive in the averaged nonlinearity. We write their sum as 


Rk(fl)=^ J2 Cki...k,kU kl ■■■V kp vi p+1 ---vZ g 6l+j‘ ..qk<5K+f...qk) , 

p<qk 1 ...kpk p+ i...kq 

where 

C( l...p N _ f 1 if OJki + . . . + CJ kp = U! k p+1 + • • • + UJkg + cj k 

< H W P+I...qkl - | o else 

This suggests to take for the effective equation the following system: 

dk = Rk(d) — 7kdk + /.ifek/3 k i k G Z d . (1-13) 


Indeed, it is proved in [7] (also see [2, 3, 4]) that, if the original equation (1.6) is well posed on time 
intervals t < 1/T, then eq. (1.13) describes the limiting behaviour of the variables a k (and, as well, the 
distribution of energy since |uk| = |uk|) in the time-scale t ~ 1/rq for any sufficiently regular initial 
data. This holds both in the presence and in the absence of the random forcing (i.e., both for /r = 0 and 
= Moreover, in the forced case we also control the limiting behaviour of the stationary solutions 
for eq. (1.6). So if the equation (1.6) and the effective equation both are mixing, then we control the 
behaviour of all solutions for (1.6) under the iterated limit lim e _>o linit-nx, . Remarkably, in this case 
the effective equation describes not only the limiting behaviour of the actions, but also that of the 
angles. I.e., it completely controls the limiting distribution of solutions. So if f(v) is a functional on 
the space of sequences v = (r>k), satisfying some mild restriction on its growth as the norm of v goes 
to infinity, and v e (t) is any solution for (1.6), then 

lim lim E/(v s (f)) -> f f(y)n{dv ), 

£—>0 t—¥ OO J 

where /x is the unique stationary measure for the effective equation (1.13) and E signifies the expectation. 
See in [3, 4, 7]. 
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1.3 Structure of resonances 

We intend to use the effective equations as a tool to investigate the energy transport in different 
physically relevant PDEs. We will show that the limiting, as e —> 0, energy transport for any specific 
equation depends on the structure of the resonances (which, in turn, is determined by the form of the 
dispersion function o>k). 

Three possibilities can occur: 

1) The resonant set is empty. Then if the degree q of the nonlinearity is even, the effective equation is 
linear. If q is odd, the equation may contain nonlinear integrable terms of the form /(/)i>k- But these 
terms do not contribute to the dynamics of the wave actions. So in any case different modes do not 
exchange energy, and no energy transport to high frequencies occurs. 

Now assume that the resonant set is not empty. We say that integer vectors ki, k 2 £ Z d are 
equivalent if there exist vectors k.s,..., k g , k £ Z d , such that the relations (1.12) hold. This equivalence 
divides Z d to clusters, formed by elements which can be joined by chains of equivalences (see [10, 11] 
for a discussion of the role of resonant clusters in weak turbulence). 

The two remaining cases are: 

2) All resonances are connected, so the whole Z d is a single cluster. In this case, in the limit when the 
volume of the space-domain goes to infinity, under some additional assumptions a new type of kinetic 
equation can be derived, the energy transport takes place and power law stationary spectra, which 
depend only on the form of the dissipation, can be obtained. 

3) All resonances are divided to non intersecting clusters. Now the energy transfer should be studied 
separately within each cluster. If sizes of the clusters are bounded, then no energy transport to high 
frequencies occurs. 

See [3, 12] for the case 1), [4, 5] for the stochastic case 2) and [13] for the deterministic case, and 
see [6] for the stochastic case 3). See [7] for discussion and for theorems, applicable in all three cases, 
deterministic and stochastic. 

Note that many examples of systems which fall to type 2) are given by equations (1.6) with com¬ 
pletely resonant spectra {tUk}, i.e. with spectra of the form cuk = a>* L?k, where I?k are integers. Av¬ 
eraging theorems for completely resonant deterministic equations (1.6) with v = 0 were discussed in 
[14, 15, 13]; also see [7]. 

Below we discuss examples for the case 2) when all resonances are connected (Section 1.4), and for 
the case 3) when the resonances make non intersecting finite clusters (Section 1.5). For more examples 
of systems of types 2) and 3) see [10]. 

1.3.1 The equations 

Our first example is the cubic NLS equation on the d-dimensional torus of size L (see [4, 5]): 

dtu — iAu = ie 2 5 |rt| 2 u + u( dissipation) + pv^random forcing) , u = u(t, x) £ C , (1.14) 

for x £ T d = R d /(27rLZ d ), where the parameter 5 = 8{L) is introduced in order to control different 
scaling for solution as the size L of the torus varies. 4 The dissipation is a linear operator, diagonal in 
the exponential base of functions on T d 

{Mx) = e iL ~ lk ' x , k£Z d }. (1.15) 

4 More exactly, (1.14) is the damped/driven cubic NLS equation. See [8, 9] for the non-perturbed NLS 
equations. 



6 S. Kuksin, A. Maiocchi 

As before, by v = (i>k,k £ Z d } we denote the the Fourier coefficients of u(x): 

u (x) = ^2 v^Mx) ■ 

k£Z d 

If d = 1, the resonance condition (1.12) takes the form 

An + &2 = &3 + fc , k\ k2 — kz k . 

All solutions for this system are such that k\ = k,k 2 = ks, or k\ = k^,k 2 = k. So in this case the 
resonant set is empty, and no energy cascade to high frequencies happens when e 2 = v — > 0. This is 
well known. 

Now consider a higher-dimensional NLS equation, write it in the Fourier variables {t>k, k £ Z d }, 
and pass to the slow time r = vt. Then, if the forcing and the dissipation are chosen in accordance 
with the prescription of the previous section (cf. ( 1 . 6 )), the equation reads 

v k = k + i5 E Uk 1 Uk 2 Wk 3 53k-7kWk + /-t&k/3 k . (1.16) 

ki ,k 2 ,k3 

Here — 7 k are the eigenvalues of the dissipation operator. The eigenvalues of the operator — A, call 
them , follow the dispersion relation 

cj k =|k| 2 /L 2 , k£Z d . (1.17) 

Below we will see that if d > 2, then the whole Z d forms a single cluster, so the equation fits the case 2). 


An interesting example of the case 3) of isolated clusters is provided by the Charney-Hasegawa- 
Mima (CHM) equation on the j5 plane (see [6] and see [8, 9] for this equation with v = 0), which we 
write as 


{—A + F)dtip — eJ(tp, Aifi) — d x il> = v( dissipation) + (random forcing) , r/) = ip(t, x) £ R . 

Here the constant F > 0 is called the Froude number and J(ip, Aip) denotes the Jacobian determinant 
of the vector (rp,Ai/j). The space-domain is a strip of horizontal size p and vertical size one, under 
double periodic boundary conditions, i.e., 

x = ( x , y) £ T 2 p l = R/ (2-KpZ) xS 1 , S 1 = M./(2nZ) . 

Again we pass to the Fourier modes 5 {uk, k = (m, n) £ Z 2 } and to the slow time r to re-write the 
equation as 


1>k = — iv 


^ifUk + , 2 2 2 1 F 2 \ E (ml + nip 2 ) 

p(m z + n z p z + F p z ) ' v 7 

ki,k2£Z 2 

x (mm 2 - mm 2 ) Vki«k 2 <5k 2 - 7k«k + ^&k/3 k > 


(1.18) 


where ki = (mi,ni), k 2 = (m. 2 ,ri 2 ) and the dispersion function has the form 


Wk = 


mp 

m 2 + n 2 p 2 + F p 2 


k = (m, n) £ Z 2 . 


(1.19) 


5 


Note that, due to the fact that the function if> is real, Vk = v! k - 
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The effective equations for eq. (1.16) and eq. (1.18) can be easily obtained on account of the general 
formula (1.13). Using it, for the NLS equation we get the effective equation 

a k = iS 52 ak 1 ak 2 Sk 3 53k<5(w JV 3k) — 7kfik + /7&k/3 k , k £ Z d , (1-20) 

ki ,k2 ,k3 

while for CHM the effective equation is the system 

5k = , , , — ! , , „ o, 52 ( m i + riip 2 ) {mm 2 - mm 2 ) 

+ (121) 
X Uk 1 Uk2 ( ^k ^(^ k ) 7k«k T /T^k/^k , k £ Z . 

It is clear that the behaviour of solutions for eqs. (1.20)-(1.21) is dictated by the structure of 
resonances since they determine the surviving terms of the nonlinearity. The geometric properties 
of the resonant set for the higher dimensional NLS equations are described in the following section, 
whereas the resonances for CHM are discussed in Section 1.3.3. 

1.3.2 Structure of resonances for the NLS equation 

In the case of 2d NLS equation each resonance is formed by four points of Z 2 which have a simple 
geometrical characterization: they form the vertices of a rectangle. Indeed, if a quadruple {k, ki, k 3 , k 2 } 
satisfies (1.12) with q = 3, then on account of the second relation we have ki — k = k 3 — k 2 . So the 
polygonal {k, ki,k 3 ,k 2 } is a parallelogram. Substituting k = ki + k 2 — k 3 in the first relation and 
using (1.17) we get 

2(k 3 • k 3 + ki • k 2 - k 2 • k 3 - ki • k 3 ) = 0 =7 (k 3 - k 2 ) • (k 3 - ki) = 0. 

I.e., k 3 — k 2 is orthogonal to k 3 — ki. So {k, ki, k 3 , k 2 } is a rectangle in Z 2 . 

It is easy to see that for any vectors k, ki £ Z 2 there is an integer rectangle of the form 
{k, ki, k 2 , k 3 }. So the equivalence, defined by the clusters of the 2d NLS equation makes Z 2 a sin¬ 
gle cluster, and the equation falls in the case 2). A graphical illustration of some resonant quadruples 
and their connections in Z 2 is displayed in fig. 1.1. 

Similar all higher-dimensional NLS equations fall in case 2). 

1.3.3 Structure of resonances for CHM 

The structure of resonances for the CHM equation depends on the shape-factor p. Below we discuss it, 
supposing for simplicity that the Froude number F is kept fixed (see [6] for the general case). We start 
with explicitly rewriting for the present case the resonance condition (1.12) (recall that ki = (mi,ni), 
k 2 = (m 2 , n 2 ) and k = (m, n))\ 

mi + 7712 = m , Til + 772 = 77 , 

mi 7772 771 (1-22) 

m\ + 7i 2 p 2 + Fp 2 + n\p 2 + Fp 2 m 2 + n 2 p 2 + Fp 2 

Solutions (ki,k 2 ,k) to these equations can be divided to different classes, according to how many 
numbers among mi, m 2 and m are non-zero: 

(i) If all three are zero, then any triad ki = (0, rii), k 2 = (0, 712 ), k = (0, ni +ti 2 ) constitutes a solution. 
As c kl ,k 2 ,k vanish in this case (see (1.21)), such triads do not form a resonance. 
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Fig. 1.1. Examples of connected resonant quadruples for the NLS equation in the Z 2 lattice. Each 
point is the vertex of at least one rectangle. 


(ii) If only one number is different from zero, then (1.22) admits no solution. 

(iii) If only one among mi, m 2 and m 3 vanishes, two subcases arise (as ki and k 2 play an exchangeable 
role): 

a. if mi = 0 (which implies m 2 = m), then = n 2 and there are two solutions, one for m = 0 , 

n 2 = n, and another for n 2 = — rii /2 = — n; 

b. if m = 0 (which implies mi = —m 2 ), then n 2 = n\ and again there are two solutions, one for 

n = 0 , ni = —n 2 , and another for ni = ri 2 = n/ 2 . 

(iv) All three are different from zero. This is the only case when the solutions depend on p. Indeed, let 
us fix a triad (ki, k 2 , k) and look for which values of p it constitutes a resonance. The second line 
of ( 1 . 22 ) may be re-written as 

00 (ki, k 2 ,k) + ai(ki, k 2 ,k, F)p 2 + a 2 (ki, k 2 , k, F)p 4 = 0 , (1-23) 

where ao,ai and a 2 are polynomials. In particular, do = m\m 2 m ( 7712 m + mim — 77117712 ). In view 
of the inequality ( x 2 + y 2 + xy) > 0, valid for nonvanishing x and y, 

ao = mim 2 m (m 2 + mj + mim 2 ) ^ 0 , 

where the use is made of the first line of (1.22). This shows that the second order polynomial in 
p 2 at the l.li.s. of (1.23) is nontrivial. So for any fixed triad (ki,k 2 ,k), where mi, m 2 and m are 
nonzero, relation (1.23) holds for at most two nonnegative values of p. 

The different types of resonances are represented in Fig. 1.2, where only the points above the 
horizontal axis are displayed (cf. footnote 5). There the resonances of type (iii) (which we will call 
standard resonances) are connected by solid lines: they form triangles symmetric with respect to the 
vertical axis m = 0, in which each point (2m, 0) is connected with (m, n) and (m, — n), for any n. The 
resonances of type (iv) (which we will call non-standard) are displayed as dashed lines: they constitute 
triangles in which none of the vertices lies on the vertical axis. 
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Fig. 1.2. Example of connected resonant triads for the CHM equation in the Z 2 lattice. Points belonging 
to different clusters of standard resonances are marked with different symbols, solid lines connect 
standard resonances, dashed lines non-standard ones. 


Since each non-standard resonance appears only for at most two values of p, then by removing (at 
most) a countable set of p’s we kill all of them. Let us denote this removed set Z C R+. The set R + \.Z 
of remaining values of p, for which no non-standard resonance appear, can be regarded as “typical”. 
Accordingly we will refer to p € R + \.Z as typical values of p (or as to the typical case). Below in 
Section 1.5 we show that in the typical case all resonances are divided to non-intersecting clusters of 
size at most 3, thus fitting the third option, considered in Section 1.3. 


1.4 NLS: the power-law energy spectrum 

Effective equation (1.20) for NLS (which, as we have seen, determines the energy spectrum) is not easy 
to handle since its completely connected resonance structure (see Section 1.3.2) makes impossible to 
split it to simpler subsystems (on contrary to the CHM equation, see Section 1.5). We present here 
a way to investigate the behaviour of solutions of (1.20) when the size L of the box goes to infinity, 
based on certain traditional for the wave turbulence heuristic approximation (see [9, 8, 16]), following 
our work [5]. This will lead us to a wave kinetic (WK) equation of the form, usually encountered in the 
wave turbulence. The treatment follows closely the paper [5], to which the reader can refer for further 
details. 

1.4.1 The limit L —>■ oo 

From the point of view of mathematics, the limit when the size L of the torus Tf, tends to infinity in 
equation (1.20) presents a serious problem. In particular, for what concerns the definition of a possible 
limiting stochastic equation. Instead of trying to resolve this difficulty, for any finite L we will study the 
expectations E(/k) of the actions Ik = ^|5k| 2 of solutions for the corresponding equation (the function 
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k h-> E(Jk) is called called the wave-action spectrum), and then pass to the limit as L —> oo only for 
these quantities. 6 

We fix L and, by making use of Ito’s formula for Ik, get from (1.20) that 


^7 k = "T 5Z - a^a^a^ak) <5 3 £<5 (w^k) 

ki ,k2 ,k3£Z^ 

— 7k/k + ^bk ^«k/3k + akflk) + bjk , k £ Z 2 . 


(1.24) 


Now we pass to the expected values, and define the moment M. 


k„ 1+1 ..:k„ 1+7l2 ( T ) o f a(r) of order m+n 2 


M , 


kl • "kn^ 
a 1 +l-" k n 1 +n 2 


(r) = (a kl • • • o kni a*k ni+1 ■ ■ ■ 5^+,^ 


(1.25) 


where (-) r stands for the expected value at time r, i.e., (f(a)) T = E(/(a(r))j for any measurable 
function /(a). Then from the system (1.24) we get 


M k k = -2 7 kA/£ + 2bl + 25 . k 6 Z d . (1.26) 

ki,k 2 ,k 3 


This system is not closed since it involves the moments of order 4. By applying again Ito’s formula, 
we can express the time derivative of moments of any order ni + n 2 as a function of the moments 
of order ni + n 2 — 2 and those of order m + n 2 + 2. The coupled system, containing the equations 
for all moments, is called the chain of moments equation (see [17]). 7 Systems of this kind are usually 
treated by approximating moments of high order by suitable functions of lower order moments in order 
to get a closed system of equations. We will show that if the quasi-stationary and quasi-Gaussian 
approximations (see below) are chosen to close the system of moment equations, then under the limit 
L —> 00 we recover a modified version of the WK equation. 

To study the sum in the r.h.s. of (1.26), we notice that if the Kronecker deltas are different from 
zero because k equals to one vector among ki, ka and k a is equal to another, then the moment is real 
and does not contribute to the sum. So we may assume that k ^ ki,k 2 , k 3 ^ ki,ka. In this case we 
calculate the fourth order moments in the r.h.s. of (1.26) through Ito’s formula (see [5]) and get 


= -(7k + 7k i + 7k 2 + 7k 3 )M k % k2 +*5 £ 

k4,k5,k6 


.*(< 


N kk 4 
J ke;k< 


,)+ 


» yk 1 k 2 k 4 rk 3 k 4 c/ JV k 3k4, a yk 2 k 5 k 6 rk 5 k 6 r , JV k 5k6 . , , 

M kk 5 k 6 OkskgOi^ k 5 k 6 J --M k k 3 k 4 OkikiOi^ k 4 k 4 J- M 1 


kr k.^kft c-k.^ki 


kk 3 k 4 


c/ JV k 5k6s\ 

kokj 5 ^ k 2 k 4 ) J 


(1.27) 


We make now the first approximation by neglecting the term containing the time derivative at the 
l.h.s. of (1.27). This can be justified, if r is large enough, by the quasi-stationary approximation (cf. 
Section 2.1.3 in [8]). Namely, let us write equation (1.27) as 


dr 


+ (7k + 7k 1 + 7k 2 + 7k 3 ) M, 


rkik 2 

J kk 3 


/• 


Notice that since all 7 k’s are positive, then the linear differential equation in the l.h.s. is exponentially 
stable. Assume that / as a function of r is almost constant during time-intervals, sufficient for relaxation 
of the differential equation. Then 

6 In the case of the non-forced equation the expectations should be taken with respect to the distribu¬ 
tion of the initial data, while for the forced equation - with respect to the distribution of the forcing 
(and, possibly, of the initial data). 

' Notice that, since the equation which we consider has a cubic nonlinearity, equations for moments 
of even order are decoupled from those for moments of odd order. 
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M, 


kik 2 

kk3 


/ 


7k + 7k i + 7k 2 + 7k 3 


We insert this in (1.26) and get 


« - 2 7k M k k + 2 bl + 25 2 V ——----- 

,, 7k + 7k i + 7k 2 + 7k 3 3 




ki,k 2 ,k 3 


Re f E + M^S^S(UJ^) 

\k4,k5,ke 


(1.28) 


We then apply the second approximation, generally accepted in the WT (see [8, 18, 9, 16]) which 
enables us to transform the previous relation to a closed equation for the second order moments. This 
consists in the quasi-Gaussian approximation, i.e., in the assumption that the higher-order moments 
(1.25) can be approximated by polynomials of the second-order moments, as if the random variables 
Vk were independent complex Gaussian variables. In particular, 


M\l$l « M\l m\ 1 M\l 


S l £ (« + <W) +-5!, 1 (<W + + % (%% + ■ (1-29) 


At this point we pass in equation (1.28), closed using the relation (1.29), to the limit L —> oo. 
This can be done by approximating sums with integrals if, instead of parametrising the modes by 
integer vectors k £ Z d , we parametrise them by vectors k = k/L from the shrunk lattice Zf, = L~ lr L d . 
Accordingly we define 



+ n 2 




l+ n 2 


7k : = 7k , h '■= fc k , 


and note that since the restriction, imposed by the Kronecker deltas, is homogeneous, then it does 
not change under this re-parametrisation. Abusing notation, we will drop the tildes in the rest of the 
Section, but will use the parametrisation by points of Z)(. 

We denote by Sk the sum, given by the second and third lines of (1.28), written in the new 
parametrisation, and note that it splits into a finite number of sums like 


st= F i( k 1 ,k 2 ,k 3 ,k 4 ,k 5 ,k 6 ). 

(k 1 ,k 2 ,k 3 ,k 4 ,k 5 ,k 6 )ez6‘ i n.E j ( 


Here is a manifold in R 6d , defined as 

E^ = {(X1,X2,X3,X4,X5,X 6 ) : Xi + x 2 = k + x 3 ,|xi| 2 + jx 2 | 2 = jkj 2 + Ix 3 1 2 

x J +X4 = x 5 + x 6 , |x J | 2 + |x 4 1 2 = I x 5 1 “ + |x 6 | 2 ,x^ = x^,x^ = x^,xf = x^j , 

where x J stands for one among the vectors k, xi, x 2 , x 3 , and {x 4 ,..., Xg } - for a permutation of the 
set {k, xi,..., X6 }\{x j }7 It is easy to see that since every F J is a regular function, then when passing 
from the sums to integrals in the limit L —> oo, each term as a function of L becomes proportional 
to L m , where m is the dimension of the manifold E^. A detailed analysis of all cases shows that the 
terms of the highest order in L in the integral correspond to terms of the form 

8 Note that the relations, defining E^ are not independent. 
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si= E 

ki,k 2 ,k 3 


ka.ks)^ 


jV k ik 2 

<H W kk 3 ) 


in the sum Sk, where k := (ki, k 2 , k 3 ) G Z 3 “ =: A4. Denote 

X k = |x = (xi,x 2 ,x 3 ) G R 3d : xi + x 2 = k + x 3 , |xi| 2 + |x 2 | 2 = |k | 2 + |x 3 | 2 | . 

This is a manifold of dimension 3d — d — 1 = 2d — 1, smooth outside the origin. The latter lies outside 
Xk if k 7 ^ 0, and is a singular point of Xk if k = 0. 

As shown in [5], in the limit L —> oo the sum Sk can be approximated by the integral 


S k 


-Pk(x) 


dx , 


J£\{ 0} < / 5 k(x) 

where <^k(x) is a certain function on Xk, smooth outside zero, such that 


(1.30) 


Vl < <Pk(x) < Vi(3d) d 1/2 , := (fi mk (mx) , _ ^ 

V2k(xi,X2,X 3 ) = V 2 k (x 2 ,Xi,X 3 ) , V?k(xi,X2,X 3 ) = lfi X3 (x 3 , X 2 , k) , 

where Vi is the volume of the 1-ball in R 2d_1 . 

By substituting (1.29) in (1.28) and using (1.30) we get the limiting (as L —> oo) equation in the 
form 


M k « - 2 7 k+ 2bi + 4d 2 L 2d ~ 1 


[ dkidk 2 dk 3 - < Pk 1 ( k ^ k 2. k 3) - 5 k^k 2 

R33\{0} 7k + 7k 1 + 7k 2 + 7k 3 3 


x <5(^1:^ 2 )(m£m£m£ + m£m£m£ - . 


Finally, we define 

n k = L d M k /2 , b k = L d/2 b k , 

(so that 53 k M k /2 —» J n k and 53k ^k —> f b k as L goes to infinity), choose 

S(L) = s 2 L 1/2 = E -—fL x/2 , 

for some e > 0, and get the kinetic equation 

h k = - 27 k n k + bk + 16e 4 [ dki dk 2 dk 3 (5 k ^ 2 (5(a; iY kk3 2 ) 

V 5 k 1 (ki, k 2 , k 3 ) / 

x -;-;-;- n kl n k2 «k 3 + n k n kl n ko - n k nk 2 % 3 - n k n kl n k3 

7k+7ki+7k 2 +7k 3 V 


(1.32) 


(1.33) 


(1.34) 


We have thus shown that, with the proper scaling of <5 and b given by (1.32)-(1.33)), the function 
n k satisfies a kinetic equation, similar to the WK equation for NLS in the classical wave turbulence 
theory (see, for instance, formula (6.81) of [9], where d = 2). The differences are two: obviously in our 
case there are forcing and dissipation, absent in the traditional WK equations. More interesting is the 
nonvanishing denominator 7k + 7k 3 +7k 2 +7k 3 which regularises the integral since it growths to infinity 
with k, and which modifies the spectra. 
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1.4.2 Power law spectra 


Now, under some additional approximation and using the well known Zakharov argument (see [8, 9, 16]) 
we will get stationary solutions of eq. (1.34) with power law energy spectra (more properly, wave-action 
spectra) {rik}. 

To do this we have to restrain our analysis to the inertial interval, i.e., to the spectral interval, 
where the damping and forcing are negligible. That is, we should consider eq. (1.34), supposing that 
the wave-vector k belongs to a sufficiently large spectral region, where the first two terms at the r.h.s. 
of (1.34) can be neglected, compared with the third. This happens, e.g., if there the solution {%} is of 
order one, while bk ^ 1 and 7 k < 1 (i.e., the damping and the dissipation are small at that spectral 
region). In the inertial interval we end up with the equation 


h k 


16 r 


dki dk 2 dk 3 5^ 2 5(w J 


^ka <ys k 1 (ki, ko, k 3 ) 


/R3d\{0} ' 7k + 7k i + 7k 2 + 7k 3 

X ^U.k 3 nk2Uk 3 T rikU-kj Uk2 UkUk2Uk3 rikllki llk 3 ) • 


(1.35) 


Notice that, while in the inertial interval we can simply approximate bk with zero, this cannot be done 
to 7k since these numbers appear in the denominator of the integral at the r.h.s. of (1.34) (their sum 
makes the denominator of the so-called collision term), and play an essential role in determining of the 
spectrum. 

The previous equation has the form of the four-wave kinetic equation (see, for instance, formula 
2.1.29 of [8]). It is well known (see [8, 9]) how to solve such equations for stationary spectra with the 
aid of the Zakharov transformations, provided that the terms 


'7“k,k 3 

'ki,k 2 


y-kHki.kz.ks) 

7k + 7k x + 7k 2 + 7k 3 


satisfy, for some ro£l, the following conditions of symmetry and homogeneity: 


T -k.ka /-j—k^.k r-i-' k,k*? rr -ki ,ko /T-Ak,Ak-? \ mr-r- k.k^ 

k!,k 2 = V,k 2 = T k 2 ,k 3 ! = ^k,k 3 , 7W,Ak 2 = A 7^4 • 

Since <p is a homogeneous function of degree 0 due to (1.31), the requirements above are met if 
on the inertial range the function 7 k can be approximated by a homogeneous function of the form 
7k = e , |k| m , where m is a real number and s' <C 1 is a small parameter to guaranty that the dissipation 
term indeed is negligible. 

We abbreviate |k| = k and look for stationary isotropic, spectra behaving as power laws of k, i.e. 
«k = nk oc k u for some real v, by searching v such that the r.h.s. of (1.35) vanishes. The result (see [5]) 
is that, in addition to the equilibrium solutions rik = C and rik = C/k 2 , which correspond, respectively, 
to the equipartition of the wave action and of the quadratic energy (Rayleigh-Jeans distribution), two 
nontrivial power law stationary distributions appear. These are the solutions: 

rik oc k K 11 , rik oc k K ” . (1.36) 


If m = 0, they coincide with the well known in the wave turbulence power-law spectra for the free NLS 
equation (1.14) 1/= o (for d, > 2, see [9]), but the dissipation modifies the power law of the decay if m ^ 0. 


1.5 CHM: resonance clustering 

Let us consider in more detail the effective equation (1.21) for the CHM equation for typical values of 
the shape-factor p. By the definition of a typical p (see Section 1.3.3), no resonances corresponding to 
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the case (iv) of Section 1.3.3 occur. We can then write the effective equation explicitly, following [ 6 ]. It 
will only involve resonances of type (iii). 

Let us consider the equations for the variables Sk with k = (m, n), separating the cases m = 0 and 
ra^O. When m — 0, the only terms which survive in the nonlinearity Rk(a) are those where ki and 
k 2 satisfy the relation (iii-b) of Section 1.3.3, while for m 7 ^ 0 only the terms falling in the case (iii-a) 
give contribution. For m = 0, the nonlinearity vanishes if n is odd, while if it is even, then 


-Rk(a) = 


p(m 2 + n 2 p 2 + F p 2 ) 


E 

m\ £Z 


2 2 
2 n p 
m 1 H-— 


®(—m i,n/2) 5 ^ 0 


which in turn vanishes because the odd symmetry in mi. On the other hand, if mi 7 ^ 0, then ki and 
k 2 are completely determined by k. So we get that 


-Rk(S) 


2 mn 


p(m 2 + n 2 p 2 + F p 2 ) 


(3 n 2 p 2 - m 2 ) Oka ( o,2n) 


(1.37) 


where we denoted k := (m, —n). Note that this formula applies for the both case m = 0 and m^O. 

Expression (1.37) entails the remarkable consequence that the hamiltonian part of the effective 
equation, i.e., the system in which forcing and dissipation are removed, 

-y-Sk = -Rk(a) , keZ 2 , (1.38) 

CLT 

is integrable and decomposes to invariant subsystems of complex dimension at most three. Indeed, if m 
or n vanish, then i?k = 0 and 5k(f) = const. Now let m,n^ 0. If 3 p 2 n 2 = m 2 , then again the equation 
for 5k trivialises. Suppose that 3 L 2 n 2 7 ^ m 2 and denote 


Ak 


2 mn 


+ n 2 p 2 + F p : 


( 0 2 2 2 \ - 
(3n p — m ) £ 


(1.39) 


Then A 5 . = —A^. Eq. (1.38) (with any fixed k) belongs to the following invariant sub-system of (1.38): 


d _ _ 

flk -^k ®(0,2n)®k? 

Clk -^k ^(0,2n)®k) 

d _ 

^“(0,2n) - 0 

(we recall that 3.(o,2n) = 0,(0 -2n) by the reality condition, see footnote 5). This system is explicitly 
soluble: if 5( 0 ,2n)(0) 7^ 0, then 


6 ( 0,2 n)(t) = Const, 

6 k (t) = 6 k (0) cos(|Ak 6 ( 0 , 2 rL)|l) + 6k(0) sgn(A k 5 ( 0 , 2 n)) sindAka^^n)^), 
where for a complex number 2 we denote 


sgn(«) = 2 /I 2 I if 2 7 ^ 0 , and sgn( 0 ) = 0 . 

The formula for a^(t) is obtained from that for 5k(f) by swapping k with k and replacing 5 (o, 2 n) by 
its complex conjugate. All these solutions are periodic, and it is easy to check that 15k | 2 + 15 ^.| 2 and 
1 6 (o, 2 n)1 2 ar e integrals of motion for eq. (1.40). 

We have established that there is no Hamiltonian exchange of energy between different modes, 
apart the coupled modes 5k and 6 k- The situation does not change much when we switch in the 
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forcing and the dissipation since the effective equation (1.21), too, splits to invariant subsystems of 
complex dimension one (if mn = 0 or 3 n 2 p 2 = m 2 ), or of dimension three (otherwise). These systems 
either are independent, or have catalytic interaction through the variables 5(o,2n), which satisfy the 
Ornstein-Ulilenbeck equation 

d - h 

^:«(0,2 n) — —J(0,2n)V(0,2n) + 1^0(0,2n)P(0,2n) ) 

and are independent from other variables. 

Being particularly interested in the exchange of energy, let us consider the equation for the actions 

I. Due to the conservation of |fik| 2 + |a|| 2 , Ito’s formula gives (cf. (1.24)) 

(^k + = — 7kAc — 7k^k + “J^k + ®k/3k) + T^k (®k/3k + «k/3k) + &k + • 

By taking the expected value, we see that the second order moments satisfy 

Mk + Ml = -7k Mi - 7k Ml + 2( 6 k + b\) , M ( ( 0 °$ = - 7 (o,„)M{°$ + 6 ? 0 , n) (1.42) 

This equations should be compared with (1.26): they show that the amount of energy contained in a 
given cluster is not transferred to other clusters and depends only on the forcing and the dissipation, 
acting in its interior. Thus the energy cascades cannot occur for typical values of p. 


1.6 Discussion 

In this Chapter we presented a method to study a weakly nonlinear PDE by investigating properties 
of the corresponding effective equation, written in terms of the nonlinearity and the resonances in the 
spectrum of the linear part of the equation. 

We have considered two examples of equations, where the structures of resonances are completely 
different. Namely, for the NLS equation all resonances are connected, and for it we have provided a 
way of getting power law stationary energy spectra (evoking a heuristic physical argument in addition 
to a rigorous mathematical theory), while for the Charney-Hasegawa-Mima equation the resonances 
form finite clusters, and for this equation we have shown (completely rigorously) that, in the typical 
case, no exchange of energy between different oscillating modes occurs. 
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